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We study the effects of losses on the entanglement created between two separate atomic gases by 
optical probing and homodyne detection of the transmitted light. The system is well-described in 
the Gaussian state formulation. Analytical results quantifying the degree of entanglement between 
the two gases are derived and compared with the entanglement in a pair of light pulses generated by 
an EPR source. For low (high) transmission losses the highest degree of entanglement is obtained 
by probing with squeezed (antisqueezed) light. In an asymmetric setup where light is only sent one 
way through the atomic samples, we find that the logarithmic negativity of entanglement attains a 
constant value — log 2 (JV) with N = 1/3 irrespectively of the loss along the transmission line. 

PACS numbers: 03.67.-a,03.67.Hk,03.67.Mn 
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INTRODUCTION 



In parallel with research on implementation and ma- 
nipulation of qubits in physical quantum systems, proto- 
cols for quantum communication and storage have been 
investigated, where the quantum information is encoded 
in the continuous quadrature variables of the electro- 
magnetic field or in the collective spin of macroscopic 
atomic ensembles. Successful teleportation entangle- 
ment and memory storage Q have been demonstrated 
making use of only Gaussian input states (coherent and 
squeezed states) and Gaussian operations (homodyne de- 
tection and bilinear interactions in the canonical con- 
jugate x and p variables of the systems). The relative 
simplicity of the implementation and theoretical analy- 
sis of Gaussian states and operations comes at a prize 
QLIEEi : distillation is provably not possible! In order 
to distil the entanglement of Gaussian states, one must 
carry out non-Gaussian operations and, at least for a 
while, leave the set of Gaussian states 6J. Although dis- 
tillation of quantum entanglement is not possible with 
Gaussian states and operations, distillation of a quan- 
tum key for cryptography has been proposed 0. Given 
the current experimental interest and the relative sim- 
plicity, it is worth while to investigate how well one may 
use Gaussian states and operations to entangle remote 
continuous variable quantum systems, coupled by a lossy 
quantum channel, and to address the related question 
about how well one may teleport an unknown quantum 
state by use of the entangled channel. 

In theory and experiment Q , entanglement between 
two trapped atomic gases is generated when the atoms 
are probed by off-resonant light (effective Faraday ro- 
tation interaction), and the transmitted probe light is 
monitored by homodyne detection. In Fig. ^b) we dis- 
play the situation we are concerned with. We imagine 
the gases to be at a large distance L from each other and 
that there is a loss of e in the light intensity between the 
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FIG. 1: Illustration of continuous variable systems sharing 
entanglement between two light pulses (panel (a)) (see, e.g., 
Q), and between two macroscopic atomic gases (panel (b)) 
(see, e.g., |2J). The losses in the transmision lines carrying 
the entanglement to large distances is indicated by the wiggly 
curves emerging from the transmission lines. 



two gases. Figure ^a) illustrates the alternative situa- 
tion where the entanglement is shared between two ligh t 
beams propagating from an EPR source 

El H US, EH 

In the present work we are primarily concerned with the 
system in panel (b), and the setup in (a) with entangled 
photons serves as a reference. 

It is an interesting observation, that the no-distillation 
theorem precludes us from establishing several such en- 
tangled ensemble pairs and concentrating the entangle- 
ment into a single system, whereas during the optical 
probing of the atomic samples, the atomic entanglement 
is indeed increasing due to the sequential transmission 
of more and more segments of the light beam, which are 
also continuous variable quantum systems. It is accord- 
ingly not clear if the entanglement achievable in a single 
pair of atomic ensembles is fundamentally limited by loss 
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and if so, what is its maximum value. We shall answer 
this question by direct computation of the state resulting 
from various measurement schemes. 

The paper is organized as follows. In Sec.^ we intro- 
duce the measure of entanglement to be used throughout 
this work. In Sec. IIIII we recall the Gaussian covariance 
matrix analysis of the continuous variable systems. In 
Sec. HVI we consider the atomic system (Fig. ^b)) in the 
asymmetric case where light is sent in the direction from 
gas T' to gas '2', only. In Sec. [V] we consider a sym- 
metrized setup with light being sent in both directions. 
In Sec. IVII we quantify the degree of entanglement that 
can be obtained in the two systems. In Sec. lVIII we con- 
clude with a discussion of the usefulness of the generated 
entanglement for teleportation with continuous quantum 
variables. 



Now, |Ai| = |A 2 | = |n - fe|, |A 2 | = |A 4 | = \n + k\ and the 
negativity N reduces to the quantity A = Var(xi — aj 2 ) — 
Var(pi + P2) known as the EPR uncertainty 0J. We 
are interested in the case of entangled symmetric states 
where 

N = A = n- k, (3) 

and N <E (0;1]. 

Since the logarithm in the evaluation of — log 2 (-/V) is a 
monotonic function it is sufficient to consider the argu- 
ment N (A) in order to quantify the degree of entangle- 
ment for two-mode Gaussian states. 



III. GAUSSIAN COVARIANCE FORMULATION 



II. MEASURES OF ENTANGLEMENT 

We are concerned with Gaussian continuous quantum 
variables systems. A Gaussian state for a vector of vari- 
ables y = (2/1,2/2, Vn) is fully characterized by its 
mean value vector m and its covariance matrix 7 with 
entrance ij given by 7^ = 2Re((j/i - - The 

degree of entanglement of a bipartite Gaussian state is 
therefore also characterized by its covariance matrix. 

In our calculations below we are concerned with a two- 
mode covariance matrix which in a local (x\,pi,X2,P2) 
basis of canonical conjugate variables reads 
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with v, 

c„, (a 



= 2Yax(xj), v Pj — 2Var(pj), and correlations 
x,p). To quantify the entanglement of such 
a state, we consider the logarithmic negativity [T2I fl3| , 
— 2k =1 log 2 [min(l, |Afc|)], which is calculated from the 
eigenvalues A^ of «ct _1 7^ , where 7^ is the partial trans- 
pose of 712 obtained by multiplying all entrances in 712 
coupling pi and other observables by —1 (i.e., c p — > — c p 
in iJTJl). and where ajk = — i[yj,Uk\ is the matrix spec- 
ifying the commutators between our x and p variables. 
Here and throughout, h = 1. For a matrix of the above 
form, we find |Ai| = IA2I and | A.3 1 = | A4 1 and only one, 
say |A| of these norms may be smaller than unity and 
hence contribute to the measure of entanglement. We 
label this quantity N — |A|, and the logarithmic negativ- 
ity is then given by — log 2 (N) . The negativity N is an 
analytical but lengthy expression of the elements of QJ, 
and we will return to special cases in Sec. V. 

In our analysis we shall also encounter the symmetric 
two-mode Gaussian state with covariance matrix 

n k 

,0 n —k I 
T12 = I 1, n „ n I ■ (2) 
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k n 
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First, as a reference, we study the influence of loss on 
the entangled fields of photons leaving an EPR-source 
(Fig. Ola)). The entanglement in these fields is max- 
imized when the entangled photons travel equal dis- 
tances L/2 so the EPR-source is placed in the center 
of the transmission lines. For this symmetric setup, the 
covariance matrix 7epr of the two beams is given by 
Eq. @ with initial values of n and k determined by the 
squeezing parameter \x\ for a pure squeezed state, i.e., 
n = cosh|x|, k — sinh|x|, and consequently, the initial 
EPR uncertainty is given by An = exp(— 2|x|)- In line 
with our previous works ]WL llq . we write Ao = 1/r, 
and refer to r as the squeezing parameter. The loss 
e' along each arm changes the canonical x variable ac- 
cording to x 1— > \/i — e'x + "v/e'iEnoisc, where x no i SC de- 
scribes the loss-induced vacuum contribution. For the p 
quadrature we likewise have p 1— > \/l — e'p + yfi p n oise- 
In terms of the covariance matrix 7epr the mapping is 
7epr ► (1 — c')7epr + f'h, where I4 is the 4x4 identity 
matrix. Therefore ri 1 — ^ (1 — e')n + e' and k 1— * (1 — e')k, 
and the EPR uncertainty, after the transmission, reads 
Aepr h- > (1 — e') An + e'. To relate to the loss e in the full 
transmission line we assume a uniform loss per propaga- 
tion distance and hence (1 — e') = \/l — e, and the final 
result for the EPR correlation is 



iEPR 



= 1 + v 7 ! - e(l/r - 1) > 1 - VT^ 



(4) 





n 



where the last expression is obtained in the limit of infi- 
nite squeezing r — > 00. The smaller the Aepr, the higher 
the degree of entanglement. 

We now turn to the case where off-resonant light 
probes two separated atomic gases. A continuous-wave 
light beam which is linearly polarized along the x di- 
rection is sent through samples of atoms with two de- 
generate Zeeman levels each described by a Pauli spin 
operator o-j (j = 1, 2, . . . , N a ) 0- As discussed in de- 
tail elsewhere 0, we describe the interaction and the 
measurement process in the time-domain. A continuous- 
wave light beam is divided into segments of duration r 
and length L = ct each of which is assumed to be short 
enough to be accurately described by a single mode of 
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FIG. 2: A beam splitter is used to model the effect of loss 
of light intensity I (1 — e)J through the transmission line 
on the (x,p) variables of the field. The output or detector 
modes (x' ,p) are determined by x — y/l — ex + y/ex v and 
p' = \/l — e P + s/epv (see text). 

the electromagnetic field. The interaction with the atoms 
and the feedback due to measurement is in turn described 
by a succession of interactions with the individual beam 
segments. For coherent light and the initially polarized 
samples, the uncertainty relations for these variables are 
minimized, implying a Gaussian state. The interaction of 
the samples with the off-resonant light beam gives rise to 
a Faraday-rotation interaction between the macroscopic 
spin operator J = l/2^j^i <T i °f the atoms and the 
Stokes operator S of the light field. 

With maximally polarized atomic samples along the x 
and — x direction, we have the classical relations J x .\ = 
J x — —J x ,2 = N a /2 where the subscripts '1' and '2' 
refer to the two gases. We assume a large number of 
photons in the probe beam such that the x component 
of the Stokes vector S x also can be treated classically. 
This leads to the introduction of the following vector of 
quantum operators 

V = (xi,Pi,x 2 ,P2,x,p) T (5) 

/ Jy,l Jz.l Jy,2 Jz,2 Sy S z ^rp 

J x \Jx J x J x \[&x \fS~x 

with canonical commutators for the three independent 
sets of modes. For a probe laser field propagating in the 
y direction and being linearly polarized in the x direction, 
the Faraday rotation interaction is oc J z ,iS z (i = 1, 2) and 
explicitly in terms of the effective Gaussian variables for 
a light segment of duration r it reads 

Hi = K T ,iPiP, (6) 

where the effective coupling constant 

■J 'N a N p h where N p h is the number of photons in each 
segment. A polarization rotation measurement on the 
optical beam, i.e., a measurement of the x variable of the 
probe light preserves the Gaussian character of the states 
and likewise does the bilinear Hamiltonian of Eq. (JBJ . 

In previous work, the effects of loss of atomic spin 
polarization and absorption of light within the atomic 
clouds were investigated 0, 0] . Here we focus on the 



loss associated with the light traveling from one gas to 
the other. This loss is partly in the transmission line and 
partly in the in-coupling devices connecting light to the 
atomic gas. The loss is modeled by a reduction (1 — e) in 
the intensity in the probe beam along the transmission 
line. We model this loss by a beam splitter which mixes 
the incoming fields x = (a + )/ \/2, p — (a — a')/ (i\/2) 
with vacuum fields x v , p v (see Fig. [SJ ■ The beam split- 
ter leads to a mixing of the operators of the incoming 
and the vacuum light modes a i— > a' — i/l — ea + i^/ea v , 
i— > a^' = \/l — eaf — iyfeal where a prime denotes 
the operators of the modes after the passage of the beam 
splitter. From these relations which apply to each polar- 
ization component of the light field, and from the defi- 
nitions of the Stokes vector components S x = a x a x /2, 
S y = (a^a_ — ala+) / (2i) , S z = (a^a+ — ala_)/2, 
the relations between the input and the primed output 
modes are easily determined. For the variables defined 
in Eq. © , we find 

and similarly for the x variables. 

The Hamiltonian describing the system when x- 
polarized light propagating in the y direction is trans- 
mitted through the two gases is H = k t ,iPiP + k-t,2P2p' ■ 
We note that the coupling strength to the second gas is 
reduced due to the loss of light intensity along the trans- 
mission line, k t ,2 = VI — e/Ct i. We denote k t = k t ,i, 
and use Eq. J7J to express the interaction Hamiltonian 
in the output modes (x' ,p' ,x' v ,p' v ) 

H = KrVl - e{p\ +P2)p - KrVePlJV ( 8 ) 

The atomic p\ and p 2 variables are coupled with equal 
strength to the output mode p', but the output vacuum 
noise term couples exclusively to p\ . 
The Gaussian variables 

V = (xi,Pi,x 2 ,P2,x / v ,p' v ,x',p') T , (9) 

describe conveniently the system with transmission loss. 
The time-evolution of the operators yi in y is described 
by Heisenberg's equations of motion, idtyi — [yi,H]. The 
evolution from t to t + t is then determined by the equa- 
tion 

y(t + r)=Sy(t), (10) 

with S fixed by the Heisenberg equation and given in 
terms of the coupling k t and the loss e. Per definition of 
the covariance matrix it follows that 

1 (t + r) = S 1 (t)S T . (11) 

Since the vacuum output components x' v ,p' v describe loss, 
we should discard them after the interaction, which sim- 
ply amounts to removing the corresponding rows and 
columns from 7, which then becomes a 6 x 6 matrix. 
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We probe the system by measuring the Faraday ro- 
tation of the probe field, i.e., by measuring the light 
field observable x' . The light field is correlated with 
the atomic samples and therefore this measurement will 
change the atomic state of the system. We denote the co- 
variance matrix before interaction and detection of light 
by 
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A C 
C T B 



(12) 



where the 4x4 matrix A describes the atomic systems, 
the 2x2 matrix B the output photon field, and the 
4x2 matrix C the correlations between the field part 
of the system that is subject to direct measurement and 
the atomic part that only feels the measurement-induced 
back-action. An instantaneous measurement of x 1 trans- 
forms a as nnium 

Aw A' = A-C{-kB-k)-C t , (13) 

where n — diag(l, 0) and where (•)"" denotes the Moore- 
Penrose pseudo-inverse. The other parts of the covari- 
ance matrix transform as B i— > 12x2, C i— > O4X2 since a 
new light segment is used in every measurement. 

In Eq. (|13|l . the change in A is proportional to k 2 which 
is proportional to the number of photons in the segment 
of duration r, and therefore proportional to r. We may 
then consider the update of the matrix A in the limit of 
infinitcsimally time steps and obtain a non-linear differ- 
ential equation with coupling constant, n 2 — k 2 /t. 



IV. ASYMMETRIC PROBING WITH 
SQUEEZED LIGHT 

Recently, we investigated the output field of an opti- 
cal parametric amplifier |2f)j | and formulated a Gaussian 
theory in the time-domain that consistently described 
the squeezing of the fields for probing times longer than 
the inverse bandwidth of squeezing. In an application 
in precision magnetometry, we showed that essentially 
the same degree of squeezing and precision could be ob- 
tained as in the ultra-broadband squeezing case where 
the squeezing is just described by the squeezing param- 
eter r. In the asymmetric case the two gases are subject 
to the interaction © , and the squeezed light is described 
by a covariance matrix diag(l/r, r), which in the output 
modes leads to Tdiag(l, 1, 1/r, r)T _1 where T describes 
the transformation from (x v ,p v ,x,p) to (x' v ,p' v , x',p') 
(see Eq. iJZJ). The atomic covariance matrix for the two 
gases is denoted by 712 as in Eq. (1). We consider the 
continuous limit of the update formulas for 712 and ob- 
tain a non- linear Ricatti differential equation (see, e.g., 
Refs. [H[ll ii The matrices in this equation 
describe squeezing of p± + p 2 , and noise introduced by 
the vacuum mode. The time-dependent covariance ma- 
trix which solves the Ricatti equation is on the form of 
Eq. QJ with 

v Xl = l + tn 2 (r + 4e(l - r) + 4e 2 (r - 1)), (14a) 



c x =t{\- e)n 2 {r + 2e(l-r)), 



tn 2 r(l - e) + (1 - r)(l - e) + r 
2tn 2 r{\ - 7) + (1 - r)(l - 7) + r'' 

tK 2 r(l - e) 
" 2tK 2 r{\ - e) + (1 - r)(l - 7) + r ' 



(14b) 



(14c) 



(14d) 



v X2 = 1 + tn 2 (l - e)(l + (1 - r)(l - e)), (14e) 



tn 2 {l - e)r +(l-r)(l-e)+r 
2tn 2 (l - Tjr + (1 - r)(l - e) + r ' 



(14f) 



We use these expressions in Sec. I VII to quantify the de- 
gree of entanglement between the two gases. We shall see 
that although some of these expressions diverge individ- 
ually for large t, the logarithmic negativity approaches a 
constant value. 



SYMMETRIC PROBING WITH SQUEEZED 
LIGHT 



We now turn to the symmetric case where equal 
amounts of noise is added to all four quadratures. The 
symmetric setup is obtained by sending light of different 
directions and polarizations through the two atomic sam- 
ples, i.e., the gases are subject to the following sequence 
of effective interactions Hi = k t \J\ — e(pi + P2)p'i — 
y/ep' vl K T pi, H 2 
II.; 

K T ^/T~~e(xi — X2)x' 4 + ^/ex' v4 k t X2, where p' iv and x' i% 
are different vacuum modes. Each of these Hamiltonians 
leads to a propagation matrix S via Heisenberg's equa- 
tions of motion. Tracing over all other degrees of free- 
dom than those of the two atomic gases and adding the 
differential equations obtained for the four Hamiltonians 
above leads to a single Ricatti equation and if we probe 
with equal strengths and equal squeezing parameter (but 
r — ► 1/r when p' components are measured instead of x') 1 
the covariance matrix is of the symmetric form J5J). The 
quantity N = A = n — k quantifies the degree of entan- 
glement, and solves the equation 



K r %/1 - e{p\ + p 2 )p' 2 - V^Pv2 K rP2, 

/t T yl — t{x\ — X2)x' 3 — y/ex' v3 K T Xi, and H± = 



d t A = a-bA 2 , 
with growth rate for the noise 

a = K 2 e(l — e + re), 



(15) 



(16) 



while the term driving the entanglement formation is 
given by 



4n 2 r(l - e) 
1 — e + re 



(17) 
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The solution to Eq. (|15fl with initial condition A(i 
0) = 1 is found by quadrature and is given by 



A(i) 



cosh(yafo<) + ■v/?sinh(va6t) 



cosh(^/abt) + j sinh(va6i) 



(18) 



In the presence of transmission line loss and squeezed 
light it follows that the steady-state value for the EPR 
uncertainty is 



A = 



1 



r(l-e) 



(19) 



Since A must be less than unity for the two gases to 
be entangled, Eq. (|19fl shows that for unsqueezed probe 
light (r = 1) it is possible to entangle the two gases only 
if e < 4/5. 

The minimum value of A (maximal entanglement) is 



and is obtained for the optimal squeezing 



1 - e 



' opt 



(20) 



(21) 



This equation has as a consequence a property that at 
first sight seems surprising: for e < 1/2 the minimum 
in A occurs for squeezed light (r > 1) as expected, 
but for e > 1/2 the minimum in A is obtained for 
antisqueezed light, r < 1. This happens because the 
squeezed light input modes contribute to the undetected 
output mode which acts as a noise term —^fep v K T px on 
the first gas in Eq. JSJ). By Heisenberg's equation of 
motion for x\, we obtain from this part of the Hamil- 
tonian the mapping x\ i— > x[ = x\ — \fep' v n T which 
means that the added noise is determined by the variance 
£K 2 Var(p^) = l/2ft 2 e(l — e + re) which is precisely half 
the factor entering the noise term in Eq. Q16[l . When we 
probe with light that is squeezed (antisqueezed) in the x 
quadrature, the p quadrature is antisqueezed (squeezed) 
and the noise term is sensitive to the noise in this anti- 
squeezed (squeezed) component 



VI. QUANTIFYING THE ENTANGLEMENT OF 
THE CHANNEL 



For the asymmetric setup with light sent only one way 
through the samples, Fig. shows a surface plot of the 
negativity N for long probing times as a function of e 
and r for r G [0.1; 10]. The negativity N is given by 
an analytical but very lengthy expression composed of 
sums of products of the variances and covariances from 

The results are summarized and compared with the 
symmetric interaction case and the EPR fields in Fig. 0] 
For a given loss, the full curve shows the optimum neg- 
ativity from Fig. [3] assuming physical limits to the avail- 
able squeezing r mul = 0.1, r max = 10. The thin full 




Loss, 8 







squeezing, r 



FIG. 3: (Color online). Negativity at asymptotic times as a 
function of loss, e, and squeezing parameter r £ [0.1; 10] for 
the atomic asymmetric case where the light is only sent one 
way through the two samples, and where only the x' quadra- 
ture is probed. 




FIG. 4: Optimal vales of the negativity for r £ [0.1; 10] at 
asymptotic times as a function of loss e through the transmis- 
sion line. The full curve corresponds to the case where the 
light is only sent one way through the two gases, and where 
only the x' quadrature is probed. The thin full curve corre- 
sponds to the case r m i n = 0.01. The long-dashed curve is for 
the symmetric case with N — e. The crossed curve is for the 
case of light emitted from an optical parametric oscillator 
and the circles denote the asymptotic form of the full curve 
in the long-time and infinite-squeezing limit r > 1, r < 1, 

ra . ran . 



curve shows the corresponding result for r min = 0.01, 
and the open circles show the results obtained with infi- 
nite squeezing: For losses less than 1/3 and in the limit 
of infinite squeezing (r > 1) we obtain 



N 



2-3e 



(22) 
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For arbitrarily high losses, we find surprisingly that the 
negativity approaches a constant for antisqueezed light 
(r < 1) ' 



N = 1/3, 



(23) 



and we find that at the particular loss of e = 1/3, 
N = 1/3 for all values of the squeezing parameter, r. 
These results mean that no matter the loss it is always 
possible (assuming that the appropriate squeezed probe 
light source is available) to entangle the two gases to a 
high level. 

The crossed curve is given by Eq. Q for the symmetric 
EPR light source operated at a finite degree of squeez- 
ing r = r max = 10. The long-dashed curve describes 
the minimum value of N for the symmetric atomic case 
(maximum entanglement) N — A m j n = e obtained by 
the optimal squeezing of Eq. I|21|) . 

From Fig.0]we see that for small loss and finite squeez- 
ing the light-atom system leads to higher entanglement 
than achievable with the EPR light source for the same 
degree of squeezing. For the asymmetric atom-light case 
we see from the full curve in Fig. 0] that the atom- 
light setup leads to a higher degree of entanglement than 
the EPR light source except in a narrow range around 
e = 1/3. In particular for losses exceeding ~ 40 %, the 
asymmetric setup presents an efficient way to generate a 
high and constant degree of entanglement. 



VII. DISCUSSION 

In summary, we have studied the effects of (light) losses 
on the achievable entanglement between two locations 
coupled by a lossy transmission line. We characterize the 
entanglement by the argument N of the logarithmic neg- 
ativity — log 2 (N) . The entanglement of atomic samples 
by optical probing is experimentally interesting because 
it works with classical optical sources and because the 
entangled states are stored in stationary media for po- 
tential later use. If only coherent light is available, the 
entanglement obtained with light being sent in both di- 
rections between the gases and with probing of suitable 
combinations of both the x and p atomic variables yields 
N = ^/e/(4 — At), and the gases become entangled as 
long as the losses are below 80 %. When squeezed light is 
available, we both have the possibility to simply prepare 
and send entangled light pulses over the transmission 
channel and we may use the squeezed light to improve 
the atomic probing schemes. For low losses, and with 
no limits to the degree of squeezing, the mere emission 
of entangled light pulses yields the largest entanglement, 
but if the degree of available optical squeezing is limited, 
it becomes advantageous to use probed atomic samples, 
and, quite surprisingly we showed that for large losses, 
a one-way probing scheme with anti-squeezed light leads 
to a finite negativity N = 1/3 and a logarithmic nega- 
tivity as large as log 2 (3) for any value of the loss. When 



the loss approaches unity, the required degree of anti- 
squeezing and probing time get larger, but as shown in 
Fig. 4., realistic squeezing levels suffice to give surpris- 
ingly small N for large losses. 

Entangled states have been proposed for various quan- 
tum communication tasks, and since losses set a prac- 
tical limit on the distance over which these tasks can 
be carried out, a scheme that attains high entanglement 
despite losses is of course interesting. To investigate if 
our entanglement is useful, we present a brief analysis of 
the achievements of the entangled gases for continuous 
variables teleportation . We have recently presented a 
general analysis of fidelities for Gaussian state transfor- 
mations, including teleportation 22], and for an entan- 
gled channel of the form (3), we reproduced the known 
fidelity for teleportation of an unknown coherent state 



F = 



1 



1 + n-k 



(24) 



This result applies to the case of symmetric probing, and 
with optimum use of squeezed probing light i|20[) the fi- 
delity gets as high as 



popt 



1 



(25) 
1/2 



approaching the classical limit of ^classical 
when the whole light field is lost. 

The entangled state obtained from the asymmetric 
probing can also be applied in the Braunstein-Kimble 
teleportation protocol, in which case the fidelity is ex- 
pressed in terms of the variances of the non-local vari- 
ables of the entangled channel p + — {p\ + pa)/\/2 and 
x- = (xi - x 2 )/\/2: 



F = 



1 



- 2Var(p< 12) ))(l + 2Var(xL 12) )) 



(26) 



The theory in Sec. IIVI provides the values Var(p^ L2 ' > ) = 

1/(2(1 + j3t)), Var(a;L 12) = (1 + act)/ 2 with = 2(1 - 
e)rn 2 / (1 + re — e) and a — n 2 e(l + re — t)/2, and instead 
of applying the protocol directly we suggest to locally 
anti-squeeze the p's, and hence p+ 2 \ and squeeze the 
x's, and hence x^ 2 \ so that V&r(p^ 2 ^) — * sVavip^} 2 ^), 



and Var(a; 



(12) N 



( 12) 

Var(x!i >)/s. The maximum value of 



(HP is obtained for s = J Var(a;L 12) ) /V'at( P ^ 2) ) , 



F 



1 -/irff i + 2 V^ 



(i->oo). (27) 



This expression is always less than or equal to 1/(1 + e), 
the optimum in the case of symmetric probing i|25|) with 
squeezed light, and equality is obtained for probing with 
squeezed/anti-squeezed light with a finite squeezing fac- 
tor, r — (1 — e)/e, which is incidentally also the optimal 
optical squeezing in the case of symmetric probing. For 
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large e, e.g., 90%, this value of r does indeed lead to 
significant entanglement, cf. Fig. 4, but the teleportation 
fidelity behaves "normally" and approaches the classical 
value of 1/2 when e — * 1. The surprisingly high en- 
tanglement does not provide a scheme for high fidelity 
teleportation! 

We note that we have not really proven that we made 
optimum use of the entangled state by the straightfor- 
ward use of the Braunstein-Kimble protocol, but we 
can give an independent argument for why (|25[1 must 
be the optimum teleportation fidelity in the case of the 
asymmetric probing scheme: Imagine a modified physi- 
cal setup where the transmission line and the second gas 
of our analysis is replaced by a loss less beam splitter 
arrangement and M separate gases, which each receive 
a fraction 1/M of the field after transmission through 
the first gas. Any one of these gases will become entan- 
gled with the first gas as described by the above theory, 
if we identify the transmitted fraction 1 — e with 1/M. 
The teleportation of an unknown coherent state is ac- 
complished by performing a joint measurement on the 
incident state and the first gas, and communicating the 
outcome to the other site, where a local operation on 
the quantum system establishes the output state. Since 
we share the entangled state with M different sites, we 
are thus able to teleport the same state to M different 
sites, i.e., produce M approximate clones of the incoming 
state. The fidelity for each of these clones, however, has 
been shown by a different argument to be limited by the 
value M/(2M - I) = 1/(1 + e) 0, which therefore must 
also be the upper limit of our teleportation fidelity. 



In conclusion, without violating fundamental results on 
the achievements of distillation, purification and cloning 
of Gaussian states, we have identified a protocol, that 
leads to finite entanglement between systems which are 
coupled by even very lossy transmission lines. If the 
probing field is split and distributed evenly between a 
number of partners, our analysis identifies instead a new 
kind of " entanglement polygamy" , where a single system 
can share finite entanglement with an unlimited number 
of other systems. Unlike a recent proposals for multi- 
partite Gaussian GHZ-type states [25|]. where the entan- 
glement can be concentrated and used for high-fidelity 
teleportation between any pair by measurements on all 
the other systems, our entanglement is already sizable, 
but not correspondingly useful for teleportation. The 
logarithmic negativity has been criticized for not always 
giving a reliable characterization of the entanglement of 
quantum systems. In our case, however, we have ob- 
served that also the more involved Gaussian entangle- 
ment of formation (2(| is finite for large losses. Rather 
than demonstrating a weakness of specific entanglement 
measures, we believe that the present study is a contribu- 
tion to our discovery of different varieties of useful and, 
possibly, useless entanglement. 
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